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ABSTRACT 

Next-generation cosmological surveys will probe ever larger volumes of the Universe, including the 
largest scales, near and beyond the horizon. On these scales, the galaxy power spectrum carries sig¬ 
natures of local primordial non-Gaussianity (PNG) and horizon-scale general relativistic (GR) effects. 
However, cosmic variance limits the detection of horizon-scale effects. Combining different surveys via 
the multi-tracer method allows us to reduce the effect down cosmic variance. This method benefits 
from large bias differences between two tracers of the underlying dark matter distribution, which sug¬ 
gests a multi-wavelength combination of large volume surveys that are planned on a similar timescale. 
We show that the combination of two contemporaneous surveys, a large neutral hydrogen intensity 
mapping survey in SKA Phase 1 and a Euclid-like photometric survey, will provide unprecedented 
constraints on PNG as well as detection of the GR effects. We forecast that the error on local PNG 
will break through the cosmic variance limit on cosmic microwave background surveys and achieve 
^■(/nl) — 1-4 — 0.5, depending on assumed priors, bias, and sky coverage. GR effects are more ro¬ 
bust to changes in the assumed fiducial model, and we forecast that they can be detected with a 
signal-to-noise of about 14. 


1. INTRODUCTION 

Upcoming cosmological surveys will probe larger vol¬ 
umes of the Universe, opening new windows to study¬ 
ing cosmological effects on horizon scales (see e.g. 
Yoo et al. 2012; Alonso et al. 2015b; Camera et al. 2015; 
Raccanelli et al. 2015). These effects include primordial 
non-Gaussianity (PNG) and general relativistic (GR) 
horizon-scale effects in the observed power spectrum. 

PNG is a key discriminator between different classes 
of inflation models. Local-type PNG (characterised by 
the parameter /nl) leaves a frozen imprint on horizon- 
scale power, allowing us to probe the primordial Uni¬ 
verse via the cosmic microwave background (CMB) and 
large-scale structure surveys. The Planck constraint 
(Planck Collaboration et al. 2015), cr(/NL) — 6.5 (using 
the large-scale structure convention), is far stronger than 
those from current galaxy surveys, but is close to the 
maximum achievable with CMB experiments, which can 
only rule out inflation models with relatively large PNG. 

Local PNG induces a scale-dependent correction to 
the bias of any dark matter tracer (Dalai et al. 2008; 
Matarrese & Verde 2008). This scale dependence can be 
probed through the two-point correlation function of the 
tracer on very large scales, allowing next-generation sur¬ 
veys to significantly improve upon the CMB constraints 
(see, e.g. Giannantonio et al. 2012; Camera et al. 2013; 
Camera et al. 2015c). 

At this level of sensitivity, neglecting GR horizon-scale 
effects would bias results. Moreover, they might hint 
at something new if GR breaks down on these scales. 
They arise via lightcone observations of dark matter trac¬ 
ers such as the number counts of galaxies (Yoo 2010; 
Challinor & Lewis 2011; Bonvin & Durrer 2011) or maps 
of intensity (e.g. the integrated 21cm signal from neu¬ 


tral hydrogen (Hi) galaxies; Hall et al. 2013), includ¬ 
ing Doppler, Sachs-Wolfe, integrated Sachs-Wolfe and 
time-delay-type terms. The lensing contribution to the 
clustering power, mediated by magnification bias, can 
also be significant on horizon scales (Alonso et al. 2015b; 
Montanari & Durrer 2015). 

Cosmic variance becomes a serious obstacle for 
horizon-scale measurements where PNG and GR sig¬ 
nals are strongest. Forecasts for next-generation surveys 
show that GR effects will not be detectable using a sin¬ 
gle tracer and PNG detection is limited to cr(/NL) > 1 
(Alonso et al. 2015b; Raccanelli et al. 2015). This calls 
for the multi-tracer technique (MT) to beat down cosmic 
variance (McDonald & Seljak 2009; Seljak 2009). 

MT has been used to explore improvements in 
the measurement of /nl (see e.g. McDonald & Seljak 
2009; Hamaus et al. 2011; Abramo & Leonard 2013; 
Ferramacho et al. 2014; Yamauchi et al. 2014). In these 
works, the lensing and GR contributions to cluster¬ 
ing power were ignored. While this may have little 
effect on o-(/nl), it can significantly bias the best-fit 
value extracted from the data (Namikawa et al. 2011; 
Camera et al. 2015b). MT has also been used to fore¬ 
cast detectability of GR effects by Yoo et al. (2012), but 
neglecting the lensing contribution and the integrated 
GR effects. Here we include all lensing and GR effects 
without making any flat-sky approximation. 

The MT technique opens a new observational window 
into probing large-scale signatures in the Universe. In 
addition to reducing cosmic variance, it also cancels the 
individual systematics of the two experiments and re¬ 
moves foreground residuals. We show here that MT is 
a game-changer in the way we design surveys to probe 
these scales, as volume is no longer the ultimate goal and 
noise reduction becomes a priority again. 
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2. THE MULTI-TRACER TECHNIQUE 

The theoretical observed fluctuations for a given dark 
matter tracer A can be written in Fourier space and New¬ 
tonian gauge in the form: 

+ /^ + ^(2^ - 


Equation (1) gives the theoretical transfer function 
Af-{z,k), while W{zi,z) is the window function cen¬ 
tered on Zi, namely the probability distribution func¬ 
tion of a source to be inside the fth bin. The prod¬ 
uct p^{z)W{zi,z) is the effective tracer’s redshift dis¬ 
tribution function inside the bin, normalised so that 
J dzp^{z)W{zi,z) = 1 for all 
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where k = (kj_,fc||) and Ti is the conformal Hub¬ 
ble parameter. The first line contains the RSD and 
leasing terms, while the next two lines constitute the 
horizon-scale GR terms. The density contrast S is 
in the comoving-synchronous gauge in order to con¬ 
sistently define the bias on large scales, and / is the 
growth rate. The correction to the Gaussian bias bo 
due to local PNG is given by Ab{z,k) = S/NLi^cC-z) — 
l]V,jnHoSc/[D{z)T{k)k^]. Here Sc — 1.69 is the critical 
matter density contrast for spherical collapse, T (fc) is the 
transfer function (normalised to 1 on large scales), and 
D (z) is the growth factor (normalised to 1 at z = 0). The 
magnification bias is Q(z, J"*) = [—cHnN^g/cHn 
where Ns{z,J- > SFA) is the background number den¬ 
sity of sources at redshift z with flux SF above the de¬ 
tection threshold The evolution bias is be{z) = 

—cHn[(l -I- z)“^A^s]/91n(l J- z). For details on the tracer- 
independent background functions E, F, G, I and J see 
Challinor & Lewis (2011). 

We see that PNG grows as k~‘^, while GR terms grow 
as 'H/k or 'H'^/k'^. The difference in the scale dependence 
and amplitude of the different terms allows GR correc¬ 
tions to be distinguished from PNG. More importantly, 
cosmic variance uncertainties come from <5, since it is a 
single realisation of the underlying probability distribu¬ 
tion. The MT technique relies on the fact that the ratio 
of different A’s is independent of cosmic variance, since 
both tracers are linear in S. 

The estimator we use is the sky map itself in the form 
of the aim- Assuming that the a£m are Gaussian, all 
the information will be encoded in the angular power 
spectrum, {airnO-^,^,) = SuiSmm'Gi. Using a Gaussian 
likelihood for the tbe corresponding Fisher matrix 
will then be enough to account for the MT effects. Ex¬ 
tending the single tracer case (Challinor & Lewis 2011) 
to multi-tracers, the angular power spectrum is given by 


Gt^ (z., ) = Jtt jdXnk Af^ (z„ fe) Af^ {zj ,k)Vcik). 

( 2 ) 

Here, Zi are the redshift bin centres and P(; is the di¬ 
mensionless power spectrum of the primordial curvature 
perturbation. The measurable transfer function in the 
bin is 

AY^{zi,k) = Jdzp^{z)W{zi,z)Af{z,k), (3) 

where p"^{z) is the redshift distribution function of tracer 
A. 


3. MULTI-WAVELENGTH SURVEYS 

In order to optimally exploit the MT method, we look 
for two surveys with significant difference in bias. We 
focus on the Hi intensity mapping (IM) survey that 
will be performed with SKA phase 1 in ‘single’ dish 
mode (Maartens et al. 2015; Santos et al. 2015) together 
with a Euclid-type galaxy survey (Laureijs et al. 2011; 
Amendola et al. 2013). We opt for the planned pho¬ 
tometric survey because it will detect a larger num¬ 
ber of galaxies than the spectroscopic option. We also 
consider variations to this survey with different noise 
and sky coverage, similar to a second-generation galaxy 
survey such as the Large Synoptic Survey Telescope 
(LSST; LSST Dark Energy Science Gollaboration 2012; 
Bacon et al. 2015). 


HI IM EXPERIMENT. In Hi IM, all galaxies with Hi con¬ 
tribute to the signal. We compute the Gaussian Hi bias, 
(z), weighting the halo bias with the Hi content in 
the dark matter halos (Santos et al. 2015). The number 
of observed sources is independent of the flux limit but 
the fractional temperature perturbation is equal to Eq. 
(1) with Q®' = 1 (Hall et al. 2013). The signal’s redshift 
distribution follows the Hi temperature, p®(z) oc Thi(z), 
which we fit using the results of Santos et al. (2015). For 
Hi IM, Afs/(l-|-z)^ = nHi(2) oc rHi(z)iJ(z)/(l-|-z)^ is the 
comoving density of Hi atoms (Hall et al. 2013), which is 
used to compute 6®. The noise angular power spectrum 
in the ith bin of frequency width Avi for an experiment 
with Nd collecting dishes, total observation time ftot and 
observed fraction of the sky /sky, is given by 


< = 


47r/skyTsys 

2NdttotAvi 


( 4 ) 


where Fys = 25 -1- 60(300 MHz/zz)^-®^ K is the system 
temperature. For SKAl, we assume Ndttot = 2 x 10® hr 
and /sky = 0.72. For Hi IM one can neglect shot noise 
(Gong et al. 2011). 


PHOTOMETRIC GALAXY SURVEY (PG). The bias, 
magnification bias, and photometric galaxies’ redshift 
distribution that we adopt are (Amendola et al. 2013; 
Raccanelli et al. 2015) 


( 5 ) 

qPG(z) = 0.2985 -b 0.5305Z - 0.1678z^ -b 0.2578z®, (6) 
pP®^(z) oc z^ exp [ - (1.412z/0.9)®/2] , (7) 


where the proportionality in the last equation is set by 
the total number of galaxies detected. For a galaxy sur¬ 
vey the noise angular power spectrum is dominated by 
shot noise, i.e.. 


■^PG “ 


-'''PG 


( 8 ) 
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Fig. 1.— Dimensionless noise and cross-noise. The solid lines 
correspond to /si;y = 0.18, dashed to /ajty = 0.36 and dotted- 
dashed to 0.72, i.e., the scenarios in Eq. (9). Note that (i) and (ii) 
have the same cross-noise since Np^ = Np^ . 

where NpQ is the number of galaxies per steradian in 
the ith bin. For the scatter between the photometric 
redshift estimate and the true redshift, we use crph( 2 :) = 
0.05(1-b 2 ) (Maet al. 2005). 

We consider three observational scenarios: (i) Euclid- 
like survey detecting 30 galaxies per square arcminute 
and covering 15,000 deg^, with 50% overlap with the 
SKAl Hi experiment; (ii) the same case but with 100% 
overlap; (Hi) a more futuristic LSST-like survey detect¬ 
ing 40 galaxies per square arcminute and covering the 
whole SKAl sky. Summarising, 

/S = (9) 

NpQ = NpQ = ^NpQ^ = SOarcmin^^. (10) 

When referring to a single tracer, we use its value of 
/sky, while in the case of MT for (i) and (ii) only the 
overlapping sky fraction is considered. Note that for (z) 
and (zz), we allow the SKA survey to be optimised for the 
smaller sky area, while maintaining the same total obser¬ 
vation time, which will make the noise decrease (Eq. (4)). 
The same cannot be done for the galaxy survey, since its 
sky coverage is assumed to be already fixed. 

HI-PG CROSS-NOISE. We also need to take into account 
the possible shot noise cross-power spectrum in each bin. 
This is due to an overlap in the halo mass range which 
the tracers probe. Even if this is small, it might be im¬ 
portant for the MT as this is the only noise showing up 
in the cross-correlation between tracers. The cross-noise 
is given by 

where dN/dM is the halo mass function, Mhi is the mass 
of Hi in a halo of mass M, and p-Hi{zi) is the Hi density. 
If the halo masses probed by the two surveys overlap, 
then 0(M) = 1; otherwise it is zero. For further details 
on the halo mass range for the SKA Hi IM survey, see 
Santos et al. (2015). The mass range for the photo -2 
survey is found by matching the number of galaxies in 


the bin given by the halo mass function with the number 
given by the redshift distribution of photometric galaxies. 

Figure 1 shows the (dimensionless) noise and cross¬ 
noise for the chosen binning. The cross-noise is only dif¬ 
ferent from zero at low redshifts when some halos will 
be detected by both surveys. While Hi IM is sensitive 
to the low halo mass galaxies, galaxy surveys are more 
sensitive to massive halos. For the bins with cross-noise, 
this is at least two orders of magnitude lower than the 
cross-angular power spectrum at the large-angle multi¬ 
poles. This might not be the case for more sensitive 
future surveys (like LSST), which may detect lower mass 
halos. 

4. FORECASTING METHODOLOGY 

We perform a Fisher analysis (Tegmark et al. 1997) for 
a set of 64 parameters 

= {/GR,/NL,lnAs,lnIlDM,ln6f,lnQj}, (12) 

with bf = bQ^Zi) (40 parameters) and Qi = (20 

parameters). We consider a bias parameter per bin, since 
due to the high precision measurements we are achieving 
with the MT technique, using a smaller number of bins 
with interpolation of parameters might impose a rather 
strong prior. 

We assume a fiducial concordance cosmology with 
Hq = 67.74km/s/Mpc, cold dark matter fraction 0 dm = 
0.26, baryon fraction Ob = 0.05, amplitude of primordial 
scalar perturbations Ag = 2.142 x 10“® and /nl = 0. 
The fiducial values for the galaxy bias and magnihcation 
bias are set by the survey specifications in § 3. The new 
parameter /gr, used as a quantiher of the measurability 
of GR effects, is defined by 

Ci = Cf+^SB+lens ^ (^^ 3 ) 

where we take /gr = 1 as the fiducial value and 
^ 5 -i-RSD-i-iens fgj. auto- and cross-correlations be¬ 

tween the density, RSD, and lensing terms. The term 
Cf^ includes all correlations that contain a GR horizon- 
scale term. We chose this definition in order to compute 
dCe/dfcB. analytically. We include all effects in Eq. (1) 
at the same time in order not to bias the accuracy on pa¬ 
rameter reconstruction (see, e.g. Namikawa et al. 2011; 
Gamera et al. 2015a,b). 

We use 20 redshift bins in the range 0 < z < 3, with 
variable size such that approximately the same num¬ 
ber of photo -2 galaxies resides in each bin. Then, we 
adopt exactly the same redshift binning for SKAl Hi 
IM, so that there is a complete overlap between the two 
tracers. We can do so thanks to the high resolution 
of an IM experiment, which allows us to tune the fre¬ 
quency windows. For the MT covariance matrix, we fol¬ 
low Ferramacho et al. (2014). We focus on large scales, 
neglecting beam effects from limited angular resolution, 
which should be negligible for small i. 

The angular power spectra for MT are computed for a 
Gaussian window function using a modified version of the 
publicly available CAMB_sources code (Ghallinor & Lewis 
2011). The code was changed to include the tracers’ 
redshift distribution and to compute 

We improve the Fisher matrix analysis to ensure nu¬ 
merical stability. Derivatives with respect to /gr and 
In bf are analytical, and numerical differentiation with 
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TABLE 1 

Marginal errors on /gr and /nl for a Euclid-like photo-2: survey. Hi intensity mapping with SKAl, and their combined MT 

ANALYSIS. 


p 

■''max 



cI/gr) 





irf/NL) 


Hi 

60 300 

PG 

60 300 

MT 

60 300 

Hi 

60 300 

PG 

60 300 

MT 

60 300 

without {ln6^,lnQi} 

1.36 

1.33 

1.58 

1.55 

0.075 

0.072 

4.57 

4.31 

5.34 

5.13 

1.23 1.12 

with {Inb/, In Qi} 

1.39 

- 

1.90 

- 

0.079 

- 

5.24 

- 

6.02 

- 

1.37 

with {In fed. In Qi} + 5% prior 

1.38 

- 

1.68 

- 

0.076 

- 

5.62 

- 

5.53 

- 

1.36 


respect to other parameters is done using a five-point 
stencil method. We use the logarithm of the parameter 
making the Fisher matrix less prone to numerical issues. 

The high dimensionality of the tomographic matrices 
and the large number of nuisance parameters (60) require 
utmost control on the matrix operations. Therefore, we 
perform matrix inversion via “inverse diagonalisation”: 
given a square matrix, its inverse is = UA“^U“^, 
where U and A are respectively the matrices of the eigen¬ 
vectors and eigenvalues of A. Thus, A is diagonal by 
construction and its inversion is trivial. This also helps 
in removing degeneracies in the Fisher matrix. Indeed, 
when marginalising over the set of nuisance parameters, 
if one or more eigenvalues (nearly) vanish, then this de¬ 
generacy does not propagate into the cosmological pa¬ 
rameters of interest (Camera et al. 2012). 

5. DISCUSSION OE RESULTS 

Unless otherwise stated, we assume configuration [i) 
for the photo-z survey, which is the more near-term sce¬ 
nario. For the Fisher matrix, we set the maximum avail¬ 
able angular scale, £min, to 2, and consider two different 
minimum angular scales, £max = 60 and 300. Table 1 
shows the Icr marginal error on /gr and /nl for the 
different tracer configurations, ^max and 60 nuisance pa¬ 
rameters. Given the high dimensionality of both the Cn 
tomographic matrix and the Fisher matrix, as well as 
the various implementations ensuring numerical stability 
of the matrix operations, the computation of the Fisher 
matrix can become unwieldy as ^max increases. There¬ 
fore in Table 1 we present results for .^max = 300 for the 
cosmological parameter set, and show the trend due to 
the inclusion of the nuisance parameters for £niax = 60 
only. The big improvement of MT over the single tracer 
is apparent: for /nl we get constraints ~ 5 times tighter, 
and for /gr the improvement is even more impressive as 
the bound shrinks by a factor > 20. Moreover, MT is 
more robust when we allow for full uncertainty on the 
bias-related nuisance parameters, as can be seen in the 
impact of a 5% prior on the nuisance parameters. 

The use of two different maximum angular multipoles 
is done not for the sake of a conservative versus opti¬ 
mistic comparison: both Umax’s are well within the lin¬ 
ear regime and the inclusion, besides Ag, of Udm and 
the nuisance parameters ensures that we do not overesti¬ 
mate the constraining power on /nl or /gr, even when 
pushing to small scales. Instead, we want to understand 
to what extent smaller scales contribute to the signal 
of PNG or GR effects, both of which are strongest on 
ultra-large scales. This also enables us to monitor the 
impact of noise. Noise usually dominates on small scales 
and is negligible in the cosmic variance limited regime, 
but Seljak (2009) and Ferramacho et al. (2014) suggested 



0.001 0.01 0.1 1 
Noise level 


Fig. 2. — Marginal Icr error on /gr (solid) and /nl (dashed) 
versus noise level for SKAl intensity mapping (blue), Euclid-like 
photo- 2 : galaxies (red) and MT (magenta) for tmax = 100. 

that the more MT is effective in ‘removing’ cosmic vari¬ 
ance, the larger the scales at which noise starts becoming 
relevant. 

In Fig. 2 we show the forecast marginal errors on /gr 
(solid) and /nl (dashed) as a function of the noise level. 
We multiply the noise of both IM and galaxy number 
counts by a fudge factor and let it vary from 0 to 1, where 
0 means a noiseless experiment and 1 is the real setting. 
As we remove noise, single tracers soon reach the cosmic 
variance limited plateau, while MT keeps improving. So, 
as cosmic variance fades, the more important the signal- 
to-noise ratio becomes. 

Seljak et al. (2009) proposed a mass-dependent weight¬ 
ing of the detected sources, which can considerably sup¬ 
press the stochasticity between halos and dark matter, 
thus reducing the shot noise contribution. By doing so, 
they showed that it will be possible for a next-generation 
Euclid-like survey to reduce the Poisson noise even by 
30%. From Fig. 2, the resulting improvement appears 
very clear (see also Hamaus et al. 2011). 

The impact of /sky and is studied by performing 
the same analysis as before, but comparing the three 
photo-z scenarios (i), (ii), and (iii) (§ 3) and changing 
from the bias of Eq. (5) to a higher bias 6 g^ = 1 + z. 
In Fig. 3 we present the corresponding forecast joint Icr 
marginal error contours in the (/gr, /nl) plane. The col¬ 
lapse from the outermost to the middle ellipses is simply 
caused by the doubling of the surveyed sky areas, but the 
innermost contours are also affected by a reduced photo¬ 
metric galaxy shot noise as the number density grows by 
30%. Also, we remind the reader that the Hi IM noise is 
linearly dependent on /sky Summarising, in Table 2 we 
quote the forecast marginal errors on the measurement 

















Hunting down horizon-scale effects with multi-wavelength surveys 


5 



Fig. 3.— Joint la marginal error contours for tmax = 300 without 
{In , In Qi} for the two photo-5: Gaussian bias models and the 
three survey scenarios. 


TABLE 2 

Marginal errors from the MT anal ysis fo r the three 

PHOTO-2: SCENARIOS WITH GAUSSIAN BIAS + Z (OR 1 + z) AND 
^max = 300. 




f^(/GR) 


f^(/NL) 

(i) 

0.071 

(0.070) 

1.12 

(0.65) 

(n) 

0.059 

(0.060) 

0.94 

(0.56) 

(m) 

0.048 

(0.053) 

0.79 

(0.48) 


of GR effects and PNG for the three photo-z scenarios 
and two biases. Results are shown for £niax = 300 and 
no nuisance parameters. 
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NOTE ADDED. While this paper was being completed, 
Alonso & Ferreira (2015) appeared, covering a similar 
topic. 
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